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Abstract: We report on numerical simulations of SU(2) lattice gauge theory with two 
flavors of light dynamical quarks in the adjoint of the gauge group. The dynamics of this 
theory is thought to be very different from QCD - the theory exhibiting conformal or near 
conformal behavior in the infrared. We make a high resolution survey of the phase diagram 
of this model in the plane of the bare coupling and quark mass on lattices of size 8^ x 16. Our 
simulations reveal a line of first order phase transitions extending from /? = to /3 = /3c ~ 2.0. 
For (5 > the phase boundary is no longer first order but continues as the locus of minimum 
meson mass. For /? > /3c we observe the pion and rho masses along the phase boundary to 
be light, independent of bare coupling and approximately degenerate. We discuss possible 
interpretations of these observations and corresponding continuum limits. 
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1. Introduction 

Non-abelian gauge theories at zero temperature and matter density can exist in a number of 
distinct phases which can be distinguished by the characteristic dependence of the potential 
energy on distance for two well separated static sources. These different behaviors of the 
potential energy can be accessed by varying the number of colors and the number of flavors 
of fermions. The collection of all of these different behaviors, when represented in the flavor- 
color space, constitutes the Phase Diagram of the given gauge theory. Up to possible dualities 
among different theories it uniquely defines each theory. In |jl| the reader will find an up to 
date review of all of the possible phases for a generic gauge theory. 

Knowing the phase diagram of strongly coupled theories has an immediate impact on the 
construction of sensible extensions of the standard model of particle interactions. Dynamical 
breaking of the electroweak symmetry is a time-honored example. It is well known that scaled 
up versions of QCD are ruled out by electroweak precision data.^ 

Using fermions in higher dimensional representations of the gauge group opens up many 
new phenomenological possibilities ^, ^. There are, in fact, a number of reasons to 
recommend using higher dimensional representations in the underlying dynamics breaking the 
electroweak theory: i) The dynamics is generally different from QCD; ii) A near conformal 
behavior can be reached for a very low number of fermions naturally reducing the contribution 

^The reader will find in ^ an exhaustive review of all of the precision data results from LEP I and II 
and how they constrain old and new models of dynamical breaking of the electroweak theory. 
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to precision observables 0]; iii) The spectrum of spin one states of these theories leads to 
interesting physical processes to be observed at the LHC 0, |, ||. 

An explicit phenomenological realization of this type of model is termed Minimal Walking 
Technicolor (MWT) Q and is based on an SU(2) gauge theory coupled to two flavors of 
adjoint quarks. This model is thought to lie close in theory space to theories with non-trivial 



infrared fixed points M, IC]. Indeed it is possible that this theory already exhibits such a fixed 



point. In the vicinity of such a zero of the beta-function the coupling constant flows slowly 
or walks. Originally such models were introduced to alleviate the flavor changing neutral 
current problem for extensions of the technicolor theory needed to give mass to the standard 
model fermions 12, |l3|, The MWT is thought to achieve such walking behavior with 
a minimal number of light (techni) quarks Q. This is the theory studied numerically in this 
paper. 

Another recent extension of the standard model which has attracted a great deal of 
interest is unparticle physics [|l5[. One simply couples a new conformal sector to the standard 
model. It is natural to identify this sector with a strongly coupled theory featuring an infrared 
fixed point. Knowledge of the phase diagram is then essential to provide natural ultraviolet 
completions of unparticle models. Making use of the analytic knowledge of the phase diagram 
one finds, for example, that it is not easy to construct gauge theories with an infrared fixed 
point able to produce spinor-type "unparticle stuff" Q. 

It is hence crucial to gain information on the phase diagram of strongly interacting gauge 
theories. Lattice methods and computational resources are now mature enough to provide 
a "first principles" systematic study of such phase diagrams, with dynamical fermions in 
the chiral limit on reasonably large lattices. Investigations of representations other than the 
fundamental have just begun |16, 17, 18, |l^, significantly extending older work on very small 
lattices |2C]. Simulation studies of many flavors in the fundamental representation have also 
become more active of late ||2l|, extending the results of |23, p4[ . 

In the current work we provide a high resolution scan of the mass/coupling phase diagram 
(not to be confused with colors/flavors/representation discussed above) of SU (2) gauge theory 
with two (Dirac) flavors of fermions in the adjoint (triplet) representation, using larger lattices 
and higher statistics than were utilized in our earlier work [^]. We find clear evidence of a 
phase boundary in the two dimensional plane of bare gauge coupling and quark mass. For 
/? < /3c ~ 2.0 the system undergoes a first order phase transition as this line is crossed. The 
latent heat of this transition goes to zero for [5 ^ (3c while the line continues to larger (3 as the 
locus of minimum meson mass. For (3 < (3c^Q see evidence for chiral symmetry breaking and 
a Goldstone behavior of the pion. Conversely, for (3 > (3c the Goldstone behavior oc niq 
disappears in a a novel way as the system is tuned close to the phase boundary, the string 
tension in lattice units is so small that we can only bound it from above, and the pion and 
rho masses drop quickly to values that are degenerate, within statistical errors. 

All this behavior is indicative of a theory that is very different from QCD, or even the 
theories with fundamental flavors that have been studied on the lattice. We will describe the 
various possible interpretations of the lattice data in the Discussion section at the end of this 
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work, but here we would like to highlight the most interesting one: all of the lattice data 
that we have obtained is consistent with either a nontrivial infrared fixed point where the 
theory becomes conformal, or a theory with an asymptotically free beta function that is very 
small for the range of scales accessible on our lattice (i.e., the coupling runs so slowly that 
we cannot start at weak coupling and still have access to the confinement scale on lattice of 
modest size). 

Either way, if this theory were to provide the mechanism of dynamical electroweak sym- 
metry breaking, the phenomenology would be radically different from a QCD-like technicolor 
scenario, and most likely naive dimensional analysis arguments would not be valid. Thus, 
the lattice results that we have obtained are quite exciting from this perspective and warrant 
further studies on larger lattices, which are currently in progress. 

In the next section we summarize some of the relevant analytical results. We then describe 
our lattice model and present our numerical results. Finally, we end with a discussion and 
interpretation of what we have found from our Monte Carlo study. A full tabulation of the 
meson masses that we have obtained is presented in Appendices A and B. 



2. Summary of the Analytical Results 

Dynamical fermion lattice simulations of higher dimensional representations are at an ex- 
ploratory stage and it is hence useful to compare the results with theoretical expectations 
obtained using various analytical methods. To gain insight one can now use, for instance, 
the conjectured all-order beta function for nonsupersymmetric gauge theories [^] together 
with the constraints from the unitarity of the conformal operators. This method constitutes 
a step forward with respect to the older approach based on the truncated Schwinger-Dyson 



equation (SD) [p5| , 2€, 27|, which is also referred to as the "ladder approximation" in the 
literature. In contrast to the ladder approximation, the all-order beta function allows one to 
determine the fermion mass anomalous dimension for any strongly coupled gauge theory at 
the infrared fixed point. Anomalous dimensions at fixed points are scheme independent since 
they represent physical quantities. The analytical phase diagram obtained by this approach, 
and a comparison of it to recent lattice results [16, [l^, 19, 22 1, is summarized in |jl|. 

In the ladder approximation the SU{2) theory with two Dirac flavors of adjoint fermions 
should be just below the conformal window where the theory develops an infrared fixed 
point Q. In the context of this approximation this means that the anomalous dimension 
of the fermion mass exceeds unity. However, according to the all-order beta function, if 
the infrared fixed point is actually reached then the anomalous dimension assumes the value 
7 = 3/4, where 7 = — dlnm/dln^u and m is the running fermion mass. If we take 7 = 1 as 
the boundary of the conformal window, the all-order beta function suggests that the SU{2) 
model is conformal in the infrared. However, the constraint coming from the unitarity allows 
7 to be as large as two before conformality is lost. Thus it is an open question whether or not 
a nontrivial infrared fixed point exists. As will be seen, the results of our lattice study suggest 
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that this fixed point may exist, though further investigations wih be required to strengthen 
the case for that conclusion. 

It is instructive to compare this theory with the case of the SU (3) gauge theory with two 
Dirac fermions in the two index symmetric representation. In "theory space" the previous 
gauge theory and the present one are very close, since the adjoint of SU (2) is equivalent to 
the two index symmetric representation. Recent lattice results [^] suggest that this theory 
may have an infrared fixed point, though more studies are needed here too. We note that 
the ladder approximation predicts that this theory is nearly conformal (i.e. walking), and 
further away from conformality then the SU{2) theory. Also, if one assumes that the theory 
is conformal in the infrared, then the all-order beta function predicts that the anomalous 
dimension of the fermion condensate is 7 = 1.3, larger than the value of 3/4 that was found 
in the SU(2) case above. If it is true that SU{3) has an infrared fixed point, it follows that 
the SU{2) theory also has an infrared fixed point, since the screening due to fermions is even 
greater in the latter case. 

As an aside, we note that it is quite interesting that for SU{2>) the anomalous dimension 
7 is larger than unity. If true, this would be quite an important result, since large anomalous 
dimensions are needed when constructing extended technicolor models that are able to account 
for the heavy quark masses, as noted in ||]. If the preliminary indications of 7 > 1 hold up 
to further scrutiny, it would overturn the common lore — but no rigorous theorem — regarding 
the anomalous dimension of the "quark" bilinear operator. 

Other interesting cases to consider are those with eight and twelve Dirac fermions in the 
fundamental representation of SU(3). The all-order beta function predicts that the conformal 
window cannot be achieved for a number of fiavors less then 8.25 (really, nine once the integer 
constraint is imposed) for the fundamental representation of SU{2>). This is confirmed by the 
latest lattice results pl| , In that work it was also suggested that the theory with twelve 
flavors has an infrared fixed point. The prediction of the anomalous dimension of the quark 
mass operator is then 7 = |. Amusingly this theory has the same anomalous dimension as the 
SU{2) two adjoint flavor theory that we study here (assuming they both possess an infrared 
fixed point). 

3. Lattice Implementation 

3.1 Action and simulation algorithm 

The lattice action we employ consists of the usual Wilson plaquette term 

= -f E E ^^Tr (u^{x)U,{x + (i)Ul{x + y)Ul{x)) , (3.1) 

with the link matrices U^{x) in the fundamental representation of SU(2), together with the 
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Wilson action for two Dirac fermions in the adjoint representation: 

= -\ E E^(^) (^/^(^) - + ^') + - U + 7^) V'lx - /u)) (3.2) 

+ J](m + 4)J]V^(x)^(x) . (3.3) 



Here adjoint links V^{x) are used, which are related to the fundamental links by 

V^\x) = ^Tr [a^U,ix)a''Ul{x)) , (3.4) 

with a"", a = 1, 2, 3 the usual Pauli matrices. 

We have simulated this theory over a range of gauge couplings /? = 1.5 — 3.0 and bare 
quark masses m ranging from —2.0 < m < 0.5 on 8^ x 16 lattices using the usual Hybrid 
Monte Carlo algorithm [28|. Typically we have generated between 400 — 2000 r = 1 HMC 



trajectories. Antiperiodic boundary conditions were used for the fermions in the time direction 
(in order to ameliorate problems with exceptional configurations at the for small quark mass) , 
whereas all other boundary conditions are periodic. 

All simuluations were run on the IBM BlueGene/L SUR machine at Rensselaer over a 
period of four months. The simulation software used is a recent, BlueGene/L architecture- 
specific version of the Columbia Physics System, modified such that SU(2) with any number 
of adjoint (Wilson or domain wall) fermions can be studied. The code has been validated by 
reproducing the results of [^] for the case of pure super- Yang- Mills. Indeed, the software was 
developed for a large-scale follow-up study of pure super- Yang- Mills that is in progress |3^ . 
The average compute rate was 70 Gflop/s, on a 128 node partition of the BlueGene/L. 

3.2 Meson operators 

We estimate the hadron masses by suitable fits to corresponding time sliced averaged corre- 
lation functions 

Goit) = E < *)roV'(x, t)i^{y, o)roV(y, o) > (3.5) 

x,y 

where Tq = 75 for the pion and Tq = 7/^, = 1, 2, 3 for the rho (the latter being averaged over 
spatial directions /i). Errors are estimated by a jackknife procedure in which fits are made 
to the meson correlators using subsets of the data, the mean and deviation of the resulting 
mass distribution yielding a mean meson mass and error. 

3.3 String tension 

We estimate the string tension as a function of lattice scale R from the large distance asymp- 
totic behavior of the Creutz ratio [BTI 



m 1 WiR,R)WiR-l,R-l) 2 

Y{R,R) = —m — ; ^ r-^ aa . 3.6 

' ^ W{R,R-l)W{R-l,R) ^ ' 
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Here W{R, R') is the expectation value of the Rx R' Wilson loop and a is the lattice spacing. 
The asymptotic behavior on the r.h.s. of ( |3.6[ ) assumes an area law for the Wilson loops. In 
practice one looks for the Creutz ratios x(-R)-R) to coalesce on an envelope where the area 
law becomes dominant. This occurs for Ra of order or larger than the scale of confinement 
£c = 1/A, where A is the usual dynamical scale of an asymptotically free gauge theory. 
In the chiral limit where the fermions are massless, the lattice spacing is a function of the 
bare lattice coupling (3 = ^/g^, through a ~ exp(— 27r^/3/6) where h is the 1-loop beta 
function coefficient. (Of course this estimate can be improved with high loop results, as 
has been considered in ^^.) Thus since a and A are physical scales, one expects to see an 
exponentially decreasing envelope for the Creutz ratios. Away from the chiral limit, there is 
a threshold mass above which the running of the coupling is altered. In that case one would 
have a = a{(3, m), where m is the bare fermion mass. Finally, we should note that since we use 
fundamental links in the Wilson loops, they are not screened by the adjoint fermions, and an 
area law emerges at scales Ra > ic, provided the theory is asymptotically free. On the other 
hand, if the theory has a nontrivial infrared fixed point, the only scale available is the finite 
extent of the lattice, x T (here, dimensionless quantities). One would therefore expect to 
see that the Creutz ratio behavior depends strongly on L, T, in contrast to what happens in 
the confining case where for aa?L'^ and aaP'T'^ very large the results become independent of 

L,r. 



4. Results 



Our results were obtained as a series of bare Wilson fermion mass scans, at fixed bare gauge 
coupling (3 = Perhaps the simplest observable to consider is the expectation value 

of the plaquette or action. Fig. |^ shows a series of scans in the "quark" mass for different 
values of (3. (In what follows we will often refer to the elementary fermions as "quarks" and 
composite states as "pions" , "rhos" , etc. Of course this is only by way of analogy, and we could 
alternatively prefix these names with "techni-".) Notice the appearance of a discontinuity for 
small /3 < /3c ~ 2.0. The data indicates that a line of first order phase transitions exists 
for small (3. Further support for this conclusion comes from the latent heat, as measured by 
the jump in the plaquette and displayed in Fig. |2|. It appears to vanish as /3 ^ /3c ~ 2.0. 
Beyond (3c we observe that the phase boundary continues as the locus of minimum pion and 
rho meson mass. The natural conclusion is that /3c ~ 2 is a second order end-point for the 
line of first order transitions. In Section |5| we will interpret the first order behavior across the 
phase boundary, at /3 < /3c, as corresponding to a bulk phase transition in the effective SU(2) 
gauge theory, in accordance with the well-known combined fundamental/adjoint plaquette 
action phase diagram |32, 33 1. 

In Fig. |3| we illustrate the behavior of the rho mass nipa by plotting it as a function 
of the bare Wilson fermion mass ma at three representative points in the phase diagram: 
j3 = 1.5, [3 = 2.0 and (3 = 2.5. (The full set of rho mass results from our studies is tabulated 
in Appendix B.) The region of quark mass to the left of the minimum corresponds to an 



- 6 - 




■ — I 


beta=1.50 


■ — I 


beta=1.75 




beta=1.90 


■ — I 


beta=1.95 




beta=2.00 




beta=2.05 




beta=2.10 


■ — ■ 


beta=2.25 


* 


beta=2.35 


■ — I 


beta=2.40 




beta=2.50 


■ — ■ 


beta=2.75 


■ — ■ 


beta=3.00 



0.5 



Figure 1: Plaquette expectation values as a function of the bare fermion mass, along lines of constant 
lattice gauge coupling P = A/g'^. It can be seen that /3c ~ 2 marks a transition, below which a first 
order phase transition is seen as the quark mass is varied. We therefore find that close to the phase 
boundary, Pc corresponds corresponds to a "bulk" transition, below which only a lattice phase exists. 
This can be understood in terms of the dynamical generation of an effective adjoint plaquette term 
in the gauge action, due to the radiative effects of nearly massless adjoint "quarks." Of course, for 
masses far enough away from the critical value the rcnormalization of the gauge action is relatively 
small and the adjoint term will not lead to a bulk transition. 



Aoki phase |34|, except that here it is for adjoint Wilson fermions, which was recently studied 
in [17p . In the case of /? > 2, as one approaches the minimum meson mass from above, there 
is a rapid drop in nip that is inconsistent with a simple linear variation with bare quark mass 
TUp oc m — TJic, as one would see in QCD, or in the case of two fundamental flavors observed 
in jl^ (cf. Fig. 5 of that reference). The data corresponds instead to a form 

mp ~ (m - mc)^/(^~'\ < e < 1. (4.1) 

Similar results are obtained for the pion, illustrated for same three values of f3 in Fig. ^ 
below. (A full tabulation of pion masses is given in Appendix A.) The figure shows the 
{mj^a)'^ as a function of the bare quark mass ma, so it is important to keep in mind that the 
dependence of the lattice spacing on /3, m also enters into the plot. There is clear evidence of 
a linear Goldstone dependence at strong coupling consistent with chiral symmetry breaking 
for P < 2.0. Conversely at P = 2.5 the data near the phase boundary line is consistent 



^Notice that the dependence of the rho mass on bare quark mass appears to depart from hnear at strong 
coupling which we attribute to the proximity of the first order phase transition 
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Figure 2: The latent heat, which appears to vanish in the /3 ^ 2 hmit. 



with a simple linear dependence of the pion mass on bare quark mass and chiral symmetry 
restoration. Again, the pion mass varies very rapidly with bare quark mass close to /3 = 2.0. 
One interpretation, which we will discuss further below, is that the lattice spacing shrinks 
significantly as one approaches the chiral limit, due to the comparable renormalization of the 
gauge coupling by gluons and and quarks in this model. I.e., the coupling walks when the 
quarks are very light, and does so over a large range of scales as the quarks approach zero 
renormalized mass. 

The behavior of the pion and rho masses along the entire phase boundary is shown in 
Fig. |5|. Two regimes are seen; a strong coupling phase with a light pion and heavy rho for 
/? < /3c ~ 2.0 and a phase for fj > fic where the pion and rho are approximately degenerate 
and independent of the bare coupling. The situation at /? ~ /3c is somewhat unclear as the 
statistical errors are large there. The phase boundary itself is shown below in Fig. ^. 

We have also made estimates of the string tension as measured by Creutz ratios x(^5 ^) 
[cf. (|3.6|) ] of various sizes. Fig ^ shows a plot of x(-R) R)-, = 1, • • • , 5, as a function of (3 as we 
move along the phase boundary rndP)- At distances of order or larger than the confinement 
scale, these ratios should coalesce on the value of aa^ where a is the string tension. For 
P = 1.9,1.95 this occurs, as it can be seen that 4) and x(5,5) coincide. For /? > 2 
the envelope where x's begin to converge cannot be seen, but the value of x(5, 5) places an 
upper bound on aa'^. It may be that much larger R values in i?) are needed, which is 
not possible on the 8^ x 16 lattice that we study here. This would be the case if the lattice 
spacing a has the very sensitive exponential dependence on /3 that would be expected from a 
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Figure 3: The "rho" mass nipa, as a function of the bare Wilson fermion mass to, for three example 
values of the bare lattice coupling f3. Note that as /3 increases past the critical value f3c « 2, the p mass 
on the phase boundary becomes small on the order of the inverse lattice size 1/L. This is consistent 
with the p becoming a massless state in the thermodynamic limit 

walking theory; for instance in the present theory using 2-loop running one would predict that 
between (3 = 2 and (3 = 2.1, aa^ would decrease by an order of magnitude and between (3 = 2 
and (3 = 2.5 it would decrease by five orders of magnitude. Given the trend in Creutz ratios 
with R aX (3 = 2, one can roughly estimate that x(7, 7) or x(8, 8) may be required before the 
envelope at this value of (3 would be seen. This would require a lattice of size 16^ X 32, which 
we are currently studying. On the other hand, and this is the possibility that we would like 
to emphasize, it could be that for (3 >2 one falls into the basin of attraction for a nontrivial 
infrared fixed point, and the area law does not hold at any scale. 

For (3 < 2 one has hints of the envelope, though the large statistical errors due to enhanced 
fluctuations at small (3 prevent us from measuring the larger loops needed for x(4, 4) and 
x(5, 5). Nevertheless, it would appear that the string tension a is of order consistent 
with a phase of the theory dominated by lattice artifacts. 

In the next section we discuss further the possible interpretations of these observations. 

5. Discussion 

Generically a lattice gauge theory will have a confining phase at strong bare coupling, and we 
believe this to be true for the present theory. Typically this is signaled by a non-zero string 
tension extracted from the asymptotic behavior of Wilson loops, or correlation functions of 
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Figure 4: The pion mass squared for example values of /3. The very sharp behavior as the bare mass 
is varied away from the phase boundary near (3 = 2 stands in contrast to the rounding that would 
normally be expected from the effects of the finite size effects. The significant decrease in the slope of 
the line as one approaches the phase boundary is presumably due to {7x1^^0)^ / {ma) ~ a, with a{(3, m) 
having a significant m dependence when the fermions are very light. This is particularly true since 
the contribution of the quarks to the running of the coupling is quite close to that of the gluons. 



Wilson/Polyakov lines. However at strong coupling the lattice theory will be dominated by 
lattice artifacts. For instance, from Fig. ^ on sees that at /3 < /3c ~ 2, the Creutz ratios 
indicate that aa^ = 0(1), so that for small /? (strong coupling) the string tension, and hence 
scale of confinement, is of the same scale as the lattice spacing a. Similarly, in Fig. ^ on sees 
that for P = 1.5 we obtain m^a > 1, indicating that the rho also lies at the ultraviolet cutoff 
scale. 

To understand whether this confining strongly coupled phase survives the continuum limit 
it is necessary to examine the behavior of, say, a Wilson loop, as the lattice spacing is sent 
to zero holding the area of the Wilson loop fixed in physical units. For a theory exhibiting 
asymptotic freedom this is accomplished by increasing p. In the case of QCD and on an 
infinite lattice this process can be continued indefinitely until we end up at the fixed point 
/3 = 00, thereby removing the ultraviolet cutoff. However it is possible that this procedure is 
interrupted by the presence of a first order phase transition at some finite bare coupling — so 
that the signal of confinement at strong coupling is not a property of the continuum theory. 
This appears to happen in this model, in the vicinity of the first order line. That is, our 
results indicate that the strong coupling phase at (3 < (3c is not continuously connected to a 
phase with a continuum limit. 
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Figure 5: Pion and rho masses along the phase boundary. Note that they become degenerate for 
(3^2. 



On the other hand, it does appear that one can move smoothly into a continuum phase 
if one starts sufficiently far away from phase boundary. In fact, this must be true since for 
large enough mass the theory is an arbitrarily good approximation of the quenched theory 
with just the Wilson plaquette action in the fundamental representation. Since that theory 
does not have a discontinuity separating the strong and weak coupling phases, we know that 
this is also true in our theory in this quenched limit. 

It is also possible that the a continuum theory with massless quarks may be obtained by 
tuning the bare quark mass in the regime (3 > Pc- The fact that the minimum meson mass 
appears to scale with the inverse lattice size is consistent with this. 

All of this can be understood in terms of radiative effects of two flavors of adjoint fermions. 
If the fermions are very light, they will generate a large adjoint plaquette term when they are 
integrated out to obtain the long wavelength effective theory. As mentioned briefly above, 
it is known that a first order transition occurs in the adjoint plaquette action theory 
and that in the mixed fundamental/adjoint plaquette theory there is a first order transition 
line when the adjoint term is sufficiently large ||3^. The interpretation of our results is 
therefore clear: if the quarks are approximately massless, the adjoint fermions lead to a large 
effective adjoint plaquette term, conventionally characterized by the coefficient (3a- As one 
passes through the phase boundary at small fundamental plaquette action coupling (3, one 
is actually moving back and forth across the first order line in the P^Pa plane. Since in the 
quenched fundamental/adjoint theory the first order line only exists for /? < 1.6, Pa > 0.7, we 



- 11 - 



-1 

-1.2 

S -1.4 
.c 

P 

-1.6 
-1.8 
-2 

1 1.5 2 2.5 3 

beta 

Figure 6: Here, the phase boundary is extracted from the minimum of rho mass at each /3. 



can understand why we too see that the first order behavior disappears for sufficiently large 
P or mass m. The fact that the transition in our theory happens instead at /3c ~ 2 would 
again be due to the radiative effects of the fermions, which will also renormalize (3. 

Since in a walking technicolor scenario we are interested in the theory with massless 
fermions, the phase of the lattice theory that is relevant for continuum physics is the phase 
where > pc- For /? — > oo we expect that the theory is driven to the asymptotically free 
fixed point known to exist in perturbation theory. However the behavior of the theory in 
the infrared is less clear. If the theory admits a new conformal fixed point then one expects 
that this will govern the long distance physics of the model and long distance features of 
the theory will be insensitive to the bare lattice coupling. In addition such a theory has 
no intrinsic scale, so that the only scales would be the lattice volume and the temperature. 
It follows that both the string tension and meson masses would scale to zero in the zero 
temperature, thermodynamic limit. This sort of behavior is certainly consistent with what 
we see in the phase (3 > (3c along the line of minumum meson mass m = m,c{P). Thus, 
our findings are consistent with the appearance of a new conformal fixed point in this theory, 
though they also leave open the possibility of a walking theory. 

However, one must be careful in drawing the conclusion that a nontrivial infrared fixed 
point exists. To take the continuum limit along the critical line requires tuning the bare 

coupling P with lattice spacing a such that finite size effects are under control. If this running 
/3(a) is sufficiently slow it can lead to extreme sensitivity in the dependence of the lattice 
spacing on bare coupling, when inverted to give a = a(/3). As was discussed in relation 
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Figure 7: Here we show the Creutz ratios along the phase boundary corresponding to minimum meson 
masses. For /3 — 1.9, 1.95 the envelope that determines aa^ can be seen from x(4, 4) and x(5j 5), since 
they coincide. For (3 > 2 the envelope where x's begin to coalesce cannot be seen, though the value 
of xi^T^) places an upper bound on aa^. For /3 < 2 there is some indication of an envelope, though 
large statistical errors prevent us from measuring the larger loops needed for x(4, 4) and x(5,5). 



to the Creutz ratio data above, small increases in the /3 would yield huge decreases in the 
lattice spacing. It is difficult to analyze such changes of scale on a relatively small lattice. In 
particular, large finite size effects can mask the true infinite volume, zero temperature physics. 
For example, the physical box size can become so small that the system deconfines and looks 
quasi-free, which would also be consistent with our data. In effect, the physics is indeed 
being dominated by a conformal fixed point — not a new infrared stable point but the usual 
infrared unstable asymptotically free fixed point. To distinguish amongst the possibilites will 
require larger lattices and a thorough study of finite size effects, so one must be cautious in 
interpreting our findings thus far. 

Another way of restating this is that any theory whose coupling runs very slowly with 
scale will necessarily generate a dynamical mass scale in lattice units (e.g., aA), which is very 
small for a weak bare coupling. To distinguish a confining theory with a small scale from 
a theory with a non-trivial infrared fixed point will then necessitate simulations on lattices 
which are significantly bigger, in lattice units, than the inverse of this small mass scale, which 
is a hard problem. And, indeed, on small lattices the physics will be governed by the usual 
ultraviolet fixed point corresponding to asymptotic freedom. On the other hand, simulations 
on larger lattices would allow us to perform a "step-scaling" analysis, in order to extrapolate 
to the infinite lattice volume behavior. We have begun studies of 16^ x 32 lattices with the 



- 13 - 



purpose of distinguishing between the walking and conformal scenarios that we have just 
described. 

Finally we would like to conclude by discussing a possible phase diagram which might be 
relevant in the situation where the theory does indeed contain a new conformal fixed point. 
Fig. m shows a cartoon of fixed points and possible RG flows for this model projected to 
the plane of bare coupling constants m). The arrows denote the flow of couplings under 
increases in length scale corresponding, for example, to a blocking transformation. The theory 
certainly contains the usual infrared unstable fixed point corresponding to (/?, ma) = (oo,0). 
A critical line corresponding to massless quarks extends out of this fixed point to smaller /? 
or stronger coupling. If a conformal fixed point exists it should form a sink for these flows as 
shown. 

In the picture we also show as a dashed line the line of first order phase transitions. Our 
data supports the conjecture that this line ends on a critical point corresponding to another 
infrared unstable fixed point. 

Furthermore, our results are consistent with the first order line and the critical line joining 
together at the critical endpoint. Any putative conformal fixed point would then serve as a 
infrared sink for massless flows out of these fixed points as shown^. 

Notice that all these fixed points are also unstable in the direction orthogonal to the 
critical line; i.e. under a mass deformation. Recognizing this fact actually allows us to draw 
a RG flow that would automatically permit a walking dynamics even in a theory inside the 
conformal window. One merely allows the theory to start near one of ultraviolet fixed points 
with a small but non-zero mass. Under blocking such a trajectory would flow initially towards 
the conformal fixed point in the vicinity of which the flow would slow before eventually flowing 
out along a direction corresponding to a mass deformation. Of course a walking scenario that 
introduces a mass for the fermions is not what is desired when trying to use the theory studied 
here for breaking the electroweak symmetry dynamically. Nevertheless, the possibility of a 
nontrivial infrared fixed point in the present theory would be interesting in its own right. 
However, there are technicolor models which make use of different gauge dynamics realizing 



this possibility as explained in [g5|, |3q| . 

In conclusion, we have found that the present theory at critical fermion mass either has a 
very slowly running gauge coupling (walking) or a nontrivial infrared fixed point (conformal). 
The behavior is drastically different from theories that do not sit near the conformal window. 
However, special difficulties emerge in the present theory due to extreme sensitivity of the 
lattice spacing on the bare gauge coupling, due to the slow running. Studies on larger lattices 
and a careful step-scaling analysis is needed, and indeed underway, in order to clarify these 
issues. 
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Appendix A: Pion masses 
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Pion masses for /? = 1.75. 



m 


pion 


error 


-1.20 


1.457 


0.011 


-1.25 


1.385 


0.012 


-1.30 


1.305 


0.010 


-1.35 


1.195 


0.011 


-1.40 


1.048 


0.015 


-1.45 


0.82 


0.42 


-1.50 


1.250 


0.026 


-1.55 


1.351 


0.016 


-1.60 


1.380 


0.0075 


-1.65 


1.386 


0.016 


-1.70 


1.360 


0.012 


-1.75 


1.3367 


0.0079 
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for (3 = 1.90. 
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Pion masses for /3 = 2.05. 
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Pion masses for /? = 2.10. 
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Pion masses for /? = 2.00. 
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Pion masses for /3 = 2.25. 
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Pion 


masses 


for /? = 2.35. 
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Pion 


masses 


for /? = 2.50. 
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for (3 = 2.40. 
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for (3 = 2.75. 
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Appendix B: Rho masses 
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for f3 = 1.50. 
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Rho masses for /? = 1.90. 
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for P = 1.75. 
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Rho masses for /? = 1.95. 
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Rho masses for /3 = 2.05. 
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Rho masses for /3 = 2.10. 
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Rho masses for f3 = 2.00. 
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Rho masses for (3 = 2.25. 
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Rho masses for /3 = 2.35. 
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Rho masses for P = 2.40. 
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0.0097 


-0.90 


1.030 


0.012 


-0.95 


0.8384 


0.0059 


-1.00 


0.605 


0.012 


-1.05 


0.3239 


0.0055 


-1.10 


0.2584 


0.0033 


-1.15 


0.5430 


0.0067 


-1.20 


0.952 


0.012 


-1.25 


1.331 


0.013 


-1.50 


1.647 


0.024 


-1.75 


1.345 


0.040 


-2.00 


0.807 


0.023 


Rho 


masses 


for P = 2.50. 



m 


rho 


error 


0.50 


2.563 


0.0057 


0.25 


2.390 


0.019 


0.00 


2.190 


0.012 


-0.25 


1.926 


0.025 


-0.50 


1.597 


0.018 


-0.65 


1.282 


0.027 


-0.70 


1.121 


0.028 


-0.75 


1.016 


0.031 


-0.80 


0.900 


0.019 


-0.85 


0.689 


0.030 


-0.90 


0.4899 


0.0078 


-0.95 


0.2935 


0.0099 


-1.00 


0.2271 


0.0073 


-1.05 


0.4207 


0.0071 


-1.10 


0.742 


0.011 


-1.15 


1.088 


0.041 


-1.20 


1.417 


0.033 


-1.25 


1.638 


0.030 


-1.30 


1.740 


0.019 


-1.35 


1.736 


0.023 


-1.50 


1.636 


0.024 


-1.75 


1.310 


0.026 


-2.00 


0.761 


0.018 



Rho masses for /3 = 2.75. 
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m 


rho 


error 


0.50 


2.490 


0.010 


0.25 


2.281 


0.011 


0.00 


2.028 


0.015 


-0.25 


1.750 


0.017 


-0.50 


1.358 


0.024 


-0.60 


1.069 


0.028 


-0.65 


0.953 


0.016 


-0.70 


0.807 


0.032 


-0.75 


0.629 


0.011 


-0.80 


0.488 


0.015 


-0.85 


0.317 


0.012 


-0.90 


0.2123 


0.0060 


-0.95 


0.305 


0.011 


-1.00 


0.527 


0.013 


-1.05 


0.835 


0.018 


-1.10 


1.147 


0.032 


-1.15 


1.465 


0.033 


-1.20 


1.652 


0.033 


-1.25 


1.738 


0.025 


-1.30 


1.764 


0.020 


-1.50 


1.621 


0.022 


-1.75 


1.301 


0.031 


-2.00 


0.681 


0.015 


Rho 


masses 


for P = 3.00. 
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